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' We study the final state interaction effects in _B — » KK decays. We find that the t channel one- 

' particle-exchange diagrams cannot enhance the branching ratios of — > K^K^ and B~ — > K^K~ 

' very sizably. For the pure annihilation process B° — > K^ K~ , the obtained branching ratio by final 

state interaction is at O(f0 
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I. INTRODUCTION 
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, B meson non-leptonic decays are important to study CP violation and to extract CKM parameters. When 

the B meson decays into two light mesons, the final state particles are energetic, so it is argued that they 
lO ' do not have enough time to get involved in soft final state interaction(FSI). In spite of the FSI, several 
factorization approaches, such as the naive factorization approach (FA) , the QCD factorization 

P^l approach (QCDF) Q, the perturbative QCD approach (PQCD) 0,01 and Soft-CoUinear-Effective-Theory 
Q^l (SCET) [2| have been established to analyze B meson decays. These approaches successfully explain many 
, phenomenons, but there are still some problems hard to explain within these frameworks, which have been 
summarized in .8j. These may be hints of the need of FSI in B decays. It has been argued that the FSI 
■ is power suppressed for the cancellation of the various intermediate states in the heavy quark limit but 
' for the finite bottom quark mass, this effect may not be very effective So FSI may be important to the 
channels that are suppressed by other factors (such as the color factor or the CKM matrix elements). For 
example, B — > KK decays are usually considered to be in the category Q| . 

FSI effects are nonperturbative in nature, so it is difficult to study in a systematic way and some different 
mechanism of the rescattering effects have been considered. In the study of D meson decays, the form factors 
are introduced to parameterize the offshellness of the exchanged particles l lUl2l , and this method still works 
in B meson case. This mechanism has been used to explain some puzzles such as -B — > tttt, ttK puzzle, 

it is argued that these puzzles can be resolved by FSI if we adopt appropriate parameters. If this is the right 
method to resolve these puzzles, it should be consistent with other channels, such as the small branching 
ratio of B ^ KK and B — > decays. The B KK decays have been measured by Belle ll^ and 
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Babar |15| , which are shown in TABLE 1( where the world average values are taken from 16]). The FA 



predictions can be consistent with the experiment for 



K°k" and S^ 



K if + if we employ the 

-rrO 



current nonperturbative inputs , thus the FSI effects may not be too large. The B ' —> K^K is a pure 
annihilation decay channel, so it is expected to be very small in FA, and the FSI can give sizable corrections. 
In this paper we will follow the method in Q , focusing on the two body intermediate states and considering 
only t-channel one-particle-exchange processes at hadron level. We will give the detailed calculation of the 
FSI effects for B KK decays in the next section, and then a brief summary in the third section. 

TABLE. 1. Measured branching fractions (xlO~^) of B KK decays 



Channel 


Babar 


Belle 


World average 


B° ^ K+K- 
B+ 'K^K+ 


1.19t°:«±0.13 

< 0.6 
1.5±0.5±0.1 


0.8 ±0.3 ±0.1 

< 0.37 
1.0± 0.4± 0.1 


U.yD_Q 24 

1.2 ±0.3 



II. FINAL STATE INTERACTIONS EFFECTS IN B ^ KK DECAYS 



Before analyzing the FSI in B KK decays, we first explore what we can get in the usual short distance 
analysis. The short distance contribution of the heavy meson decays can be expressed in terms of some types 
of quark diagrams: V , the penguin emission diagram; £ , W-exchange diagram; A, W^-annihilation diagram; 
Va, the penguin annihilation diagram (space- like); Vew, the electroweak penguin diagram; V, the vertical 
W loop diagram (time-like penguin). The penguin dominated B KK decays can be expressed as: 



A{B- 
A{^ 



K°K^') 
> K^K-] 
K+K-] 



V- 

= V 
= £ + V 



Va - -^Vew - 
- Va — -^Vew 



V, 



A, 



(1) 



In factorization approach, there is no emission tree diagram contribution to these decays. The annihilation 
diagrams A,£,V,Va are power suppressed which can be neglected in the calculation. They are usually 
believed to be long distance dominant. So the short distance amplitudes read: 



A{b" - K°K") = ,,^f^F,^^{ml){ml - ml)[VubV:,(a-, + r^al) + V.,V:Aal + r^al) 



(2) 



and A{B- K°K-) = A(b" K°K°) , ^(s" ^ K+K-) = 0, where Vub,Vud,Vcb and Vcd are CKM 
matrix elments, = 2'm%/[mb{ms + mq)]. a^''^ are combination of Wilson coefficients for four quark 
operators defined in in ref.j^. 

ai ^ Ci + ■^C'i+i) (* = odd) 



ai ^ Ci + gC**-!! (* = even) 



(3) 
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From quark-hadron duality, the decay amplitude can be got from either quark picture or hadron picture. 
The result should be equal. However, neither of the two pictures are fully understood in the B decays. 
The factorization theorem tells us to calculate the short distance contribution perturbatively and the long 
distance parts using hadronic picture. Thus a double counting problem may arise. To avoid double counting, 
we adopt leading order Wilson coefficient at the scale for naive factorization approach instead of QCDF 
(which includes some virtual corrections from long distance) for short distance calculations of i? ^ KK. 

When we calculate the long distance contributions to the decays, we consider only the CKM most fa- 
vored two body intermediate states, such as D^*^ D^*\'k-k, pp. The quark level B ^^{pp) KK di- 
agrams are shown in Figure 1. We can see that the this diagram has the same topology as the penguin 
diagram or W^-exchange diagram. From Eq.(l), we can see that this kind of diagrams can contribute to 
B K°'K",K+K-,K°K- simuhaneously. When the intermediate state is d'-*^+ D'^*^- {D'-*^+ D'-*^"), only 
penguin topology works, so it cannot contribute to the K^K" decay. 



d d u 




d d d 



FIG. 1: Quark level diagram for B tt+tt- K'^KO{K+K-) 

The hadron level diagrams are given in Figure 2. We focus on the t channel one-particle-exchange processes, 
furthermore, we consider only the case that the two intermediate particles are on shell, i.e. we only keep the 
absorptive part of diagrams in Figure 2, which gives the main contribution. 




c d 



FIG. 2: Hadron level diagrams for long distance t channel contribution to -B — > KK 

The absorptive part of the diagrams in Figure 2 can be calculated with the following formula: 

X T*{p3P4-^piP2), (4) 



which can be deduced using the optical theorem 

Taken FSI corrections into account, the topological amplitudes are: 



V = VsD + iAbs{a + b + c + d), 

£=iAbs{a + b). (5) 



Then the decay amplitudes turn to: 



^ K^K^) =V + Vew + iAbs{a + b + c + d), 
A{B- K°K-) = V + Vew + iAbs{a + 6 + c + d), 

A{B^ K+K-) = iAbs{a + b). (6) 



To perform the calculation, we introduce the relevant Lagrangian density 



0: 



Ci = --^Tr[F^,{V)F'^''{V)] + ^gvppTr{V^'P 5^ P) + gvvpt''''''^Tr{d^,V,d^VpP), (7) 

Cd = -igD'DpiD'd^^P^.Df - D*;d^'P,jD^^) - \gD'D'pe^.o.0Dl^'d'' P'^ 9^" of^ 

-igDDvD] D\vn] - 2fD*Dve^..Md''Vn]{Dl IT D*^^ - of^ D^) 
MgD'D'vD*"^ D,{vn} + ^ifD-D-vDllid^^V - O'^V^Dy , (8) 

where P and are pseudoscalar and vector multiplets respectively. Here we take the convention e^^^"^ = 1. 

Using Eq. and the Feynman rules derived from the Eqs. Q and ijHJl, we can get the leading long 
distance rescattering amplitude: 

Abs{a)= / '-^ 9l'KnA{,B -^Ti+TT )- S^^i' 

with 

A(B"^^+7r-) = *^A<-(m2)(m|-m2)[KbK"<i(«i+< + «io+^xK + «8)) 
+ ycbV:a{al + a'i^ + rl{al + al))], 

H, = -bl ■ P2 + P3 ■ P4 + PI ■ P4 + P2 ■ P3) - ' "''l^"'^ " , (10) 

where we denote the momentum by B{pb) Tr{pi)'!r{p2) K(jp^)K{p4), 6 is the angle between pi and ps, 
and r!^ — 2m'^/[mb{mu + md)]- Here F{t,mK*) is the form factor introduced to denote ofFshellness of the 
exchanged particle, which is usually parameterized as |s|: 

F{t,m)^{^;^-^r. (11) 

It is normalized to unity at i = rri^ ( t is the invariant mass of the exchanged particle), where we usually take 
n = 1. The cutoff A should not be far from the physical mass of the exchanged particle, where we choose 

A = mexc + ri^QCD- (12) 
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The parameter r/ depends not only on exchanged particle, but also on the external particles involved in the 
strong interaction. If it is determined from the B — > tttt branching ratios, then we can employ it in B — > KK 
decays for SU (3) symmetry. 

Likewise, the absorptive parts of the other diagrams are given by 



Abs{b{K)) 

Abs{b{K*)) 

Abs{c) 
Absid{Ds)) 

Abs{d{D:)) 



|pi|(icos0 



m 



K 



r 2A^''{m?) 



.Gf 



IQ-KUlB 



t — m%. 



^ |pi|dcos0 2 F'^it^fnoJ 



IGirniB 



QDsD'K' 



t — mi 



.Gf 



{mB+mD')A^^'{mjj,)H5 



2Ai'^\ml.)_^, 



rriB + mo* 



n 

I IQwrriB 



2Ai^'{ml,)^, 



mB + mo* 



(13) 
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where 



rr I \ "Pi- PSPI - PA P2- P3P2 'PA , Pi' P2P1 ' P3P2 ' PA 

H2 = (Pa ■ Pa) 2 2 ^ ' 

TO J TO2 TO^TOj 

TT/ r N (Pi -Ps)!??! •Pb)(P4 -Pi?) {P2 ■ Pa){P2 ■ Pb){p3 ■ Pb) 

H2 = P3 ■ PS)(P4 • pb) 2 2 

TOf TO2 

(Pi ■Pz){p2 ■Pa){pi ■Pb){p2 - pb) 



H3 = 2{pi-P4){p2-P3)-2{pi-P2)(p3-Pa), 



H3 = m%[{pi ■ Pi){p2 ■ Ps) - {Pi ■ P2){P3 ■ Pa)] + {pi ■ Pb){p2 ■ Pb){p3 ■ Pa) 

-{P2 ■Pb){p3 ■Pb){pi - pa) - {pi ■Pb){pa ■Pb){p2 ■ P3) + {P3 ■Pb){pa ■Pb){pi ■ P2), 

Ha, = -(p3 -pi) + 



{Pi ■P3 - ml){ml~p2 - Pi) 



H5 = {P3 ■ Pa) - 



(Pi ■P3){pi - Pa) _ {P2 ■P3){P2 - Pa) {Pi ■ P3){P2 ■Pa){pi ■ P2) 

9 ? 2 9' 

toJ TOj TO^TOj 



jTi ( u X (Pl ■ PB)iPA ■ Pb){P1 ■ P3) 

H5 ^ [P3 ■ Pb)[pa ■ pb) 2 

TOf 

(P2 ■Pb){P3 ■Pb){p2 - Pa) {pi ■Pb){p2 ■ Pb){pi ■ P3){P2 - Pa) 

9 ' 99 ' 

TOJ m|m2 

He = 2{pi ■ P2){P3 ■ Pa) - 2{pi ■ P4){p2 ■ Ps), 

H'e = ni%[{pi-p2){p3-PA)-{pi-pA){p2-p3)]-{Pl-PB){p2-pB){p3-PA) 

+ {P2 ■Pb){p3 ■Pb){pi - pa) + {pi ■Pb){pa ■Pb){p2 ■ Ps) - {P3 ■Pb){pa ■ Pb){pi ■P2)- (14) 



and 



^ D+D-) = t^foF^^^'imDiml ~ ml)[V,,V:,{a, +al + a^, + to,,/tob« + a^)) 

+ VubV:Ml+<^ + mD/mB{a'^,+al))\ (15) 

To proceed the numerical calculation, we use the parameters as follows: the Fermi constant Gp — 
1.166 X IQ-^GeV-^; the CKM matrix elements Vcb = 0.041, = -0.224, |Kb| = 0.0037, Kd = 0.974; The 
phase angle 7 — 60°; the meson and quark masses niB ~ 5.279GeV,mK ~ 0.4986*6^, m;, — AAGeVjiris = 
OmGeV.md = 0.004Gel^; the decay constants U = 0.132Gey, fo = 0.20GeV, fp = 0.2l6GeV, 
Jd" = 0.23Gey, U< = 0.16Gey; The form factors are from the light-front model [isj: i^^^^(O) = 0.35, 
Af''(0) = 0.22, Af''(0) = 0.20, F^^{m%) = 0.68, A^^' (771%,) = 0.65. The coupling relevant to the 
K*Kt: can be extracted from the K* Kn experiments: gK*+K°■K+—^■^^ a-nd we take gpKK — 4.28 
and gpKK- — 8\/2 0]. The coupling of D*DK and DID*K can be related to go'Drc by SU{3) symme- 
try. In this work we neglect the SU(3)symmetry breaking effect and employ the coupling as gwDK = 
y/rniymD^ g D* D' K — SD'Ott = 17.9. Similarly, we also use the symmetry to determine the parameter rj in 
the form factor, where the best fit from the B ttK decay is rjp — ^^(.j^^c*)-) = 0.69 0], in this work we 
choose r/ = (0.8, 1.0, 1.2) x 0.69 to include the SU{3) breaking effect. 

The rescattering effects can produce the strong phases, it may change the CP asymmetry behavior of short 
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distance calculation. The time dependent CP asymmetry of B° K^K^ is defined as 

T{B {t) ROR ) + T{BO{t) ^ R^TC) 
= Aj^ok" cos(AMt) + Sj^oK° sin(AMi), 



(16) 



with AM the mass difference of the two mass eigenstates of neutral mesons. And the direct CP asymmetry 
and the mixing induced CP asymmetry parameters are defined as, 

2Im(Aj^(,Ko) 



+ 1' 



(17) 



where the corresponding factor A^o^" ~ e~^''^^. 

Using the theoretical inputs mentioned above, we get flavor-averaged branching ratios for the short distance 
contribution as 



B{B° R^lC) = 0.94 X 10-*^, 
B{B+ lfR+) = 1.0 X 10"^ 



(18) 



And there is no direct CP violation since there is only one kind of contribution (pure penguin). After consid- 
ering rescattering effects, things will change, since more contributions with different phases are introduced. 
We summarize our numerical results in TABLE 2. 

TABLE 2. CP averaged branching ratios and CP asymmetries of S — > RR decays 



Channel 


r?(x0.69) 


Branching ratio(x 10 ^) 


Akk 


Skk 




0.8 


0.99 


-0.03 


-0.03 


BO ^ ROjf 


1.0 


1.1 


-0.04 


-0.04 




1.2 


1.2 


-0.06 


-0.05 




0.8 


0.009 


-0.04 


-0.56 


BO R+R- 


1.0 


0.021 


-0.04 


-0.55 




1.2 


0.042 


-0.03 


-0.55 




0.8 


1.1 


0.10 




B+ ^ R°R+ 


1.0 


1.2 


0.14 






1.2 


1.3 


0.18 





Prom this table, we can see that the FSI cannot enhance the branching ratio of B^{B ) —^ R^R sizably 
because the FSI increase (decrease) the real part for B^ — > R^T^t^B^ — > R^r'^)^ but decrease(increase) the 
imaginary part. The total effects don't make the average branching ratio change much. As the parameter 
T] gets larger, the FSI effects become more important and the larger strong phase is produced, so the 
absolute value of direct and the mixing induced asymmetry increases. For the charged B meson decays, 
the FSI effects are more important for Figure 2 (a, b) give double contribution (due to the interchange of 
the intermediate particles). So contrary to B^ — > R^R^ case, the direct CP asymmetry becomes positive. 
The B^^B^) — > R+R~ results are purely from the FSI effects, its branching ratio are of the order O{10~^), 



which is consistent with PQCD prediction [lOj in quark diagram calculation. It seems to be a proof for quark 
hadron duaUty. The D{D*)D{D*) intermediate states cannot contribute to B^{B^) K+K^ through t 
channel processes, the strong phase of this channel comes from the Wilson coefficients, so the calculation 
gives a small direct CP asymmetry. 

In ref 8], the DD tttt annihilation diagrams which have the same topology with vertical W loop 
diagrams, are introduced to resolve B tttt puzzle. It gives an dispersive part which can reduce B'^ — )■ tt+tt^ 
branching ratio as well as enhance B'^ n'^n'^ one. Considering SU{3) symmetry, these diagrams can 
contribute to B ^ KK at the same level as i? ^ tttt, we quote their results here (in units of GeV): 
DisA = 1.5 X 10^^VcbV*j — 6.7 x 10^''VubV*jj. If we consider this effect in B ^ KK case, the branching 
ratio for B K^K^ is enhanced to about 2 x 10^^, while the B'^ K'^k'^ branching ratio is reduced to 
about 6 X 10~^, which is not favored by B ^ KK experimental data. 

The B KK decays have also been calculated with the QCD factorization^^ and PQCD approach Q], 
in which part of the long-distance effects has been included. These methods depend strongly on theoretical 
inputs, such as the chiral factor(or equivalently, the current quark mass), so they also give large error. The 
QCDF calculations give (branching ratios are CP averaged, also for H2()|l ^: 

r/rO j:fO-Tf>\ _ -, or;+0.41+0. 70+0. 13+1. 09 , r)-6 
0{n 1\ J\ ) — -L-O0_o. 36-0. 48-0. 15-0. 45 ^ ' 

n(fl- K^K-^ - 1 Qfi+0. 45+0. 72+0. 14+0.91 -, r,-6 
K(R" K+K-\ — n ni-i+O-OOS+0. 008+0. 000+0.087 -,r)-6 

oyn j\ j\ ) — L).uio_g 005-0.005-0.000-0.011 ^ ' 

Acp{B- ^ K'K-) = -16.3t^:™tl™ X 10-^ (19) 
AcpiB" ^ K^k") = -l6-7tt-'7ttit\-'rtii x IQ-^. 



And the PQCD calculations give: 



B{B" K^lC) = 1.75 X 10-^ 
B{B- -> K°K-) = 1.66 X 10"^ 
B{B" K+K-) = 0.046 X 10"^, 
AcpiB- K"K-) = 0.11, 
Acp(S° i^°X°) = 0, 
Acp{B° K+K-) = 0.29. 



(20) 



For the branching ratio, with the error, all the calculations can be consistent. As for the CP asymmetry, 
PQCD and QCDF have opposite sign, our calculation is consistent with PQCD for B^ K^K^ , while our 
results have the same sign with QCDF for B^ K^k'^ . More experimental data are needed to test these 
predictions. 
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III. SUMMARY 

In this paper we study the FSI effects in B — > KK decays. We find that if we consider only the dominant 
t channel one-particle-exchange diagrams, the FSI effects cannot change the branching ratio of 
and B^{B~) ^ K^{K^K~) sizably, which is consistent with the current experimental data. We also 
predict the branching ratio of the B^{B^) — > K^K~ at O(10~^) by purely t channel FSI, which is consistent 
with the PQCD prediction. We also calculate the CP asymmetry in the B — > KK decays. We test the DD 
annihilation diagram (which is of great importance to resolve B ^ ww puzzle in FSI) contribution and find 
it not favored by B — > KK data. 
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